Single photons from Pb+Pb collisions at CERN SPS, QGP vs. hadronic gas 
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In a hydrodynamic model, we have analyzed the direct 
photon data obtained by the WA98 collaboration in 158 A 
GeV Pb+Pb collisions at CERN SPS. The transverse expan- 
sion of the system was taken into account. Two scenarios, 
(i) formation of quark-gluon plasma and (ii) formation of hot 
hadronic gas, were considered. Both the scenarios describe 
the data equally well. However, hadronic gas scenario require 
very high initial temperature (~ 300 MeV) and it is difficult 
to conceive existence of hadron gas at that high temperature. 
If the hadronic fluid has small radial velocity (0.2c-0.3c) ini- 
tially, the data are well explained in the hadronic gas scenario 
with reasonable initial temperatures. 

PACS numbers(s):12.38.Mh,13.85.Qk,24.85.+p,25.75.-q 



Recently WA98 collaboration has published their sin- 
gle photon emission data for 158 A GeV Pb+Pb collisions 
at CERN SPS @. Much interest was aroused after the 
publication of the WA80 preliminary results of the S+Au 
single photon data Q , as it was hoped that the they can 
be a conclusive probe of the much debated quark-gluon 
plasma (QGP), expected to be produced in relativistic 
heavy ion collisions. The preliminary data were analyzed 
by several authors. Xiong and Shuryak || analyzed the 
data assuming a mixed phase formation and found ex- 
cess photons. Srivastava and Sinha Q analyzed the data 
considering two possible scenarios after the collision, one 
with the phase transition to QGP, the other without it. 
It was claimed that the data were explained only in the 
phase transition scenario. We had also analyzed thepre- 
liminary versions of the WA80 direct photon data Q . It 
was shown that formation of viscous hadron gas in the 
initial state, can explain the data. The revised version 
of the data H were also analyzed by several authors in- 
cluding us 0. It was concluded that the data arc not 
sensitive enough to discriminate between the two alter- 
nate pictures, e.g. formation of quark-gluon plasma and 
formation of hot hadronic gas. 

Several authors have analyzed the recent WA98 sin- 
gle photon data [gj-|l~i|| . Peressounko et al JllJ concluded 
that the data could be explained only with a small initial 
radial velocity. Also the data could not distinguish be- 
tween a QGP and a hot hadronic gas scenario. Srivastava 
et al (|] could fit the data in the QGP scenari, without 
initial radial velocity. But the two loop photon rate used 
by them was not corrected for a factor of 4. They did 
not considered the hadronic scenario. They argued that 
in the pure hadronic scenario, initial temperature of the 
hadronic fluid will be large. Hadronic density will be ~ 
10 hadrons/J m 3 . It is unphysical to consider hadronic 



gas at such a high density. Initial Hadronic scenario was 
considered in Q. It was concluded that only hadronic 
model with medium modification (hadrons were formed 
with zero mass) could explain the data. 

In the present paper we analyze the WA98 single pho- 
ton data in the no phase transition (NPT) scenario. A 
hot hadronic gas is assumed to be formed in the initial 
state. It expands, cools and freezes out at freeze-out tem- 
perature (If). It will be shown that WA98 single pho- 
ton data could be well explained in this scenario, with 
reasonable hadron density, if one assumes a small ini- 
tial fluid velocity. In fill initial velocity distribution was 
assumed to be linearly increasing. We assume that at 
the beginning of one fluid hydrodynamic stage, initial 
radial velocity follows the distribution of initial energy 
density, which we assume to be of Woods-Saxon form, 
with radius and diffuseness parameters, R = 6.4/rn and 
a = 0.54/ m. In [[ll] hadronic equation of state comprises 
all the hadrons in the particle data book. We choose 
to use hadrons up to mass 2 GeV. The cut off is arbi- 
trary. The resulting equation of sate is reasonably well 
described by pu — <x/jT 4 , with a h — 59.57r 2 /90. Fig.l 
compares the analytic expression with numerical results. 
Analytic expression overestimate the pressure in the tem- 
perature region of 200-400 MeV. However, we have used 
the analytic expression only to obtained the initial tem- 
perature of hadron gas, and the error introduced due to 
the approximation is around 10-15%. To be complete, we 
also analyze the data in the phase transition (PT) sce- 
nario, when QGP is formed in the initial state. The equa- 
tion of state for QGP was assumed to be p q = a q T 4 — B 
with a q = 42.25tt 2 /90. The bag constant B was obtained 
from the Gibbs condition p QGP (T c ) = phad(T c ). 

We solve the hydrodynamic equations d^T^ — 
in 3+1 dimension assuming cylindrical symmetry and 
boost-invariance in the longitudinal direction. The rele- 
vant equations are well known jl^,|l3| and are not repro- 
duced here. The inputs of the hydrodynamic equations 
are the initial energy density or temperature (Tj) and ra- 
dial velocity (v™) at (proper) time Tj. Tj is the thermal- 
isation time beyond which hydrodynamics became appli- 
cable. For a given t% the initial temperature Tj of the 
fluid (hadronic gas or QGP) can be obtained by relating 
the entropy density with the observed pion multiplicity 
(assuming pion decoupling to be adiabatic) |Q, 
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where c = 27r 4 /45C(3) and Ra is the transverse radius 
of the system (assumed to be 6.4 fm for Pb+Pb colli- 



sions). 6 = corresponds to central collisions. In table 
1., we have shown the initial temperatures as obtained 
from eq.j] in the two considered scenarios, the QGP and 
the hot hadronic gas. dn/dY was assumed to be 750 [||. 
It can be seen, in both the scenarios, initial temperatures 
are comparable. Corresponding hadron densities are also 
shown in table 1. For n=0.2 and 0.4 fm, it is ~ 10 and 2.9 
hadrons per fm~ 3 respectively. It is unlikely than at such 
high density hadrons can retain their identity. For larger 
Tj, initial temperatures are comparatively small and the 
densities have acceptable values. The other parameter 
for the hydrodynamic evolution with transverse expan- 
sion is the initial radial velocity v™ % . It is customary to 
assume that the initial v* m =0. However, it is possible 
that the fluid (QGP or hadronic gas) possess some small 
radial velocity at initial time Tj. As will be shown here, 
initial small radial velocity can affect the photon spectra 
considerably. 

For the single photons from hadronic gas we include 
the following processes, 

(a) 7T7T — > (Tf, (b) irp — > 7T7, (c) Ll> — > 7T7, (d) p — > 7T7T7 
(c) Tip — > Al — > 7T7 

rates for which are well known |l^Ji^ j. 

Rate of production of hard photons from QGP were 
evaluated by Kapusta et al [O. To one loop order, 



E 



dR 



I 



d?p 2tt 2 



aa s ^)T 2 e- E/T H- CE 
f 



y a s T 



(2) 



where the constant c ~ 0.23. The summation runs over 
the flavors of the quarks and e/ is the electric charge of 
the quarks in units of char ge o f the electron. 

Recently Aurenche et al evaluated the production 
of photons in a QGP. At two loops level Bremsstrahlung 
photons {qq(g) — > qq{g)l found to be dominating the 
compton and annihilation photons. However, their cal- 
culations overestimated the photon yield by a factor of 4 

The rate of production of 

photons due to Bremsstrahlung (corrected for the factor 
of 4) was evaluated by them as, 
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where Jt ~ 4.45 and Jl ~ —4.26 for two flavors and 
3 colors of quarks. For 3 flavor quarks, Jt ~ 4.8 and 
J L 4.52. I{E,T) stands for, 
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and the poly-logarith functions Li are given by, 



Aurenche et al |18| also calculated the contribution of 
the qq with scattering, which was also overestimated by 
the same factor of 4. The corrected rate is, 
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We would like to mention that two loop photon rate 
from QGP is not complete. Higher loops contribute to 
the same order |2^|. Also Landau-Migdal-Pomeranchuk 
effect been neglected p^l ■ 

Before we present our results we would like to make a 
brief comment on the direct QCD photons. Direct QCD 
photons are produced from the early hard collisions of 
partons in the nuclei and in Pb+Pb collisions make sig- 
nificant contribution to the high pt yield. Gallmeister et 
al p0| claimed that prompt photons are able to explain 
the high pt data in Pb+Pb collisions. Dumitru et al |23j | 
also arrived at a similar conclusion including the nuclear 
broadening effects. However, this point is still controver- 
sial due to uncertainties in prompt photon emission at 
AA collisions. Thus Alam et al |J and also Srivastava 
and Sinha Q calculated the prompt photon emission for 
Pb+Pb collisions. It was seen that for pr >2 GeV, direct 
QCD photons alone can describe the data within a factor 
3-8 only. 

We first present the photon spectra obtained in the 
phase transition scenario. As in ref. M we assume the 
critical temperature to be T C =I80 MeV. Data were found 
to be insensitive to the exact value of T c B. The ini- 
tial radial velocity (v™ 1 ) was assumed to be zero. In 
fig.2, computed photon spectra for different initial times 
t,;=.2,.4,.6 and .8 fm, with one+two loop order rates are 
shown. In the inset of Fig.l, we have compared the pho- 
ton yield obtained with one loop rate and one+two loop 
rate. For initial time Tj = 0.2 fm, two loop contribution 
is more than 50% at large pr- Higher loop contributing 
to the same order, photon rates, are not complete. But 
for Tj > .2/m, two loop contribution is quite less and the 
rates can be regarded more or less complete. We find 
that with corrected photon rate, QGP scenario, unlike 
in ref. ||], do not describe the data well. In ||, pho- 
ton rates were not corrected for the factor of 4. With 
correct photon rate, yield is reduced by a factor of 2, 
compared to uncorrected rate. Consequently, while Sri- 
vastava et el H found very good fit to data for Tj-0.2 fm, 
presently data are underestimated by the similar factor of 
2. For higher thermahsation times, the data are further 
underestimated particularly at high px side. The initial 
temperatures are low enough to produce requisite num- 
ber of high pt photons. In fig.2, the dotted line shows 
the contribution of direct QCD photons, as calculated by 
Alam et al ||. For Tj > 0.2 fm, the hard QCD photons 
contribute more than the thermal photons. It is evident 
that for thermahsation time T;=0.2 fm, thermal photons 
together with hard QCD photons describe the data satis- 
factorily. For higher thermahsation times, thermal pho- 
tons together with hard QCD photons underpredict the 



data. We do not elaborate on this scenario. Just a few 
comments are in order. The equation of state of the 
hadronic sector now consists of hadrons with mass less 
than 2 GeV, while in ref. || hadrons with mass less than 
2.5 GeV were included. The other difference is the ini- 
tial energy density profile. Srivastava and Sinha assumed 
the initial energy density profile to follow the wounded- 
nucleon distribution, while we have used the standard 
Woods-Saxon profile. Despite these differences, our re- 
sults are very similar to the results obtained in || . The 
data are not sensitive enough to the details of the calcu- 
lations. The data can not distinguish whether hadronic 
sector comprised with hadron with mass less than 2.5 
GeV or with mass less than 2 GeV. Also, the data are 
insensitive to the details of initial energy density profile 
so long they are not very different. May be at RH1C or 
LHC energy, data will be sensitive on these details. 

The results obtained in the no phase transition sce- 
nario, when hadronic gas is assumed to be formed in the 
initial state are presented in fig. 3. We have presented 
photon spectra for different initial times, Ti=.2,.4,.6,.8 
fm. The results are similar to those obtained in the 
QGP scenario. It is not surprising. As told earlier, 
with resonance hadronic gas, the degrees of freedom are 
comparable to the QGP. Thus initial temperatures are 
nearly same either in hadron gas or in QGP. It was 
noted quite early that hard photon production rate are 
nearly same in hadronic gas or in QGP 0. Here also 
we have shown the contribution of hard QCD photons 
(dotted line). As in the scenario with QGP formation, 
here also for Ti > 0.2/m hard QCD photons contribute 
more to high px photon yield than the thermal photons. 
For Ti=0.2 fm, hadronic gas scenario describe the data 
very well. For higher thermalisation times, the descrip- 
tion gets poorer especially in the high px sector. Hard 
QCD photons fill up substantially in the high px sector, 
but as stated earlier they are less by a factor of 3-8. It 
would seem that the WA98 single photon data equally 
well described in PT as well as NPT scenario. However, 
existence of hadronic gas at a temperature of 300 MeV is 
extremely unlikely. At this temperature, hadron density 
is very large ~ 10/m~ 3 . It is unlikely that at such a high 
density hadrons can retain their identity. Hadronic gas 
scenario at such a high density is unphysical. For higher 
thermalisation time 7";=0.6 fm, when the initial tempera- 
ture and density are 211 MeV and 1.3 hadrons per fm~ 3 , 
though the scenario is physical, the data gets underpre- 
dicted by a factor of 3-8 (taking into consideration of 
hard QCD photons). One of the deficiency of the model 
is the neglect of pre-equilibrium photons, emitted before 
Ti. Traxler and Thoma p4| calculated pre-equilibrium 
photons and found them order of magnitude less than the 
equilibrium photons. Roy et al J|5| also calculated the 
pre-equilibrium photons. They used Fokker-Plank equa- 
tions and found that pre-equilibrium photons are less or 
at best equal to the equilibrium photons. The results 
refers to RHIC and LHC energies. At SPS one expects 
still lower contribution of pre-equilibrium photons. At 



SPS, the system is far from chemical equilibrium and 
number of quarks and anti-quarks will be less. Even if 
pre-equilibrium photons contribute as much as equilib- 
rium photons, the WA98 data will remain underpredicted 
for thermalisation time Ti ~ 0.6 fm or more. 

It would seem that though in NPT scenario good de- 
scription of data is obtained, physical consideration (i.e. 
very high hadron density) will render that picture un- 
acceptable. In the calculation presented till now, the 
initial radial velocity (u*™ 4 ) of the fluid was assumed to 
be zero. However it is possible that at initial time Tj the 
fluid has some small radial velocity v l ™ . Several authors 
|p[-|rT| have used initial velocity to fit the direct photon 
data. Source of v 1 / 11 may be the collisions among the 
constituents, which lead to the local equilibrium. Peres- 
sounko et al |ll| used the high Px component above 2 
GeV in the ir° spectrum to argue for it. However, argu- 
ment in favor of initial fluid velocity is weak. High px 
part of the spectrum might have come from hard pro- 
cesses pq| . Also several hydrodynamical calculations, 
with out initial velocity, could explain a host of experi- 
mental data p^-|29]| . Apparently large px data require a 
initial fluid velocity. 

In fig. 4, we have shown the photon spectra in 
the hadronic gas scenario with initial fluid velocity 
t>™=0,.l,.2,.3 (in units of c) for initial temperature of 
Ti=211 MeV corresponding to initial time Ti=0.6 fm. 
vl m had considerable effect on photon spectra. It en- 
hances the px- Good fit to data is obtained for v™ 1 
=0.3c. It is also obvious that with v'" u in the ranges 
of 0.2 - 0.3c, it will be possible to fit the WA98 data in 
a hadronic gas scenario with physically acceptable initial 
time and temperature. 

Good fit to the data with small initial velocity bring 
back the hadronic gas scenario into contention. It is no 
longer possible to say that the WA98 data indicate quark- 
gluon plasma formation only. 

Initial fluid velocity will also affect the photon spectra 
in the phase transition scenario. As shown here in this 
scenario, WA98 data are underpredicted with initial time 
Ti > 0.2 fm. It will be possible to fit the data with r, >0.2 
fm, if small initial velocity is assumed. In ref. U goof fit 
to data was obtained for r^=l fm and initial velocity in 
the ranges 0.2-0. 3c. 

To summarize, we have analyzed the recent WA98 
single photon data using a hydrodynamic model. Two 
scenarios were considered, the phase transition scenario 
where a QGP is formed in the initial state, and the 
no phase transition scenario where hot hadronic gas is 
formed initially. Both the scenarios gave good descrip- 
tion to the data. QGP scenario require that the initial 
time and temperature of the QGP is 0.2 fm and 340 MeV 
respectively. The hadronic gas scenario also require an 
initial time of 0.2 fm and temperature of 304 MeV. As the 
hadron density is very large at this temperature, it would 
seem that the data is described by QGP only. However, 
it was shown that with small initial radial velocity in the 
range 0.2c-0.3c, the WA98 data can be well described 



in the hadronic scenario with reasonable initial tempera- 
ture. The present analysis thus suggests that the WA98 
single photon data are not conclusive. It can not dis- 
criminate between two alternate scenarios currently in 
vogue. 
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TABLE I. The initial temperature of the QGP and the hot 
hadronic gas for different initial times t;. Also shown are the 
corresponding hadron density 



n (fm) 


T? GP (MeV) 


T? ad (MeV) 


P? ad (fm- 3 ) 


0.2 


341 


304 


10.41 


0.4 


271 


242 


2.89 


0.6 


237 


211 


1.32 


0.8 


215 


192 


0.78 


1.0 


200 


178 


0.52 



70 
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FIG. 1. Pressure as a function of temperature for the 
hadronic gas comprising hadrons with mass less than 2 GeV. 
The solid line is a fit using a^T 4 , an = 59.5. 




12 3 4 



Pj (GeV) 

FIG. 3. The single photon yield in the no phase transi- 
tion scenario for four different n's, .2, .4, .6, .8 fm (from top 
to bottom). Corresponding tempertures are listed in table 1. 
Experimental points are also shown. The dotted line is the 
direct QCD photons calculated in ref. I| 
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FIG. 2. The single photon yield in the phase transition 
scenario for four different initial times, Tj's, .2, .4, .6, .8 fm (from 
top to bottom). Corresponding temperatures are listed in 
table 1. Experimental points are also shown. The dotted line 
is the direct QCD photons calculated in ref. In the inset, 
we have compared the photon yield obtained with one loop 
rate (dashed lines) with the one+two loop rate. 
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FIG. 4. The single photon yield in the no phase transition 
scenario for four different initial radial velocityu^ nl =0,.l,.2,3 
(in units of c) . The initial time and temperatures are Ti=.6 fm 
and Ti=211 MeV respectively. The dotted line is the direct 
QCD photons calculated in ref. ||. 



